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We show that the mechanical properties of a worm-like-chain (WLC) polymer, of contour length L 
and persistence length A such that t = L/\ ~ 0(1), depend both on the ensemble and the constraint 
on end-orientations. In the Helmholtz ensemble, multiple minima in the free energy near t = 4 
persists for all kinds of orientational boundary conditions. The qualitative features of projected 
probability distribution of end to end vector depend crucially on the embedding dimensions. A 
mapping of the WLC model, to a quantum particle moving on the surface of an unit sphere, is 
used to obtain the statistical and mechanical properties of the polymer under various boundary 
conditions and ensembles. The results show excellent agreement with Monte-Carlo simulations. 
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I. INTRODUCTION 

Microtubules and actin polymers constitute the struc- 
ture of cytoskeleton that gives shape, strength and motil- 
ity to most of the living cells. They are semiflexible poly- 
mers in the sense that their persistence lengths A are of 
the order of their chain lengths (the statistical contour 
lengths) L such that the stiffness parameter t = L/X is 
small and finite. For example, actin, microtubule and 
double stranded DNA (dsDNA) have A = 16.7 Atw[ll,[2, 
5.2 mmp] and 50 nrnfSi] respectively. In physiological sit- 
uation L of dsDNA inside a cell may vary in between mil- 
limeters to a meter with average length in human being 
~ 5cm, whereas typical contour lengths of microtubules 
can be ~ 10/Lt?7i[3]. The contour length of actin filament 
can be as large as 100Aim[l|. In the in vitro experiments, 
the contour lengths of bio-polymers can be tailored chem- 
ically; e.g. in the experiment described in Ref.[l| the con- 
tour lengths of actin polymer have a distribution up to 
L — 30/im. For a polyelectrolyte like DNA the persis- 
tence length A can also be tuned a little by changing the 
salt concentration of the medium. The relevant parame- 
ter in deciding the mechanical properties is the stiffness 
parameter t, contour length measured in units of persis- 
tence length. While it is obvious that in the thermody- 
namic limit of i — !■ cx), the Gibbs (constant force) and the 
Helmholtz (constant extension) ensemble predict identi- 
cal properties, the same is not true for real semiflexible 
polymers which are far away from this limit. In biological 
cells actin filaments remain dispersed throughout the cy- 
toplasm with higher concentration in the cortex region, 
just beneath the plasma membrane, microtubules, on 
the other hand, have one end attached to a microtubule- 
organizing centre, centrosome, in animal cells. Thus bio- 
logically important polymers may float freely or may have 
one of their ends fixed. Even the end orientations of poly- 
mers play a crucial role in many important phenomena. 
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For instance, microtubule-associated proteins attach one 
or both of their ends to microtubules to arrange them 
in microtubule bundles [4]. Again, in gene-regulation 
often DNA-binding proteins loop DNA with fixed end 
orientations 5, 6, 7]. Thus it becomes important to un- 
derstand the statistics and the mechanical properties of 
semiflexible polymers with different possibilities of end 
orientations and ensembles. 

During the last decade many single molecule experi- 
ments have been performed on semiflexible polymers 0, 
H, These have been done by using the optical 

tweezers 0, the magnetic tweezers and the AFMsfT^. 
In the optical tweezer experiments one end of a polymer 
is attached to a dielectric bead which is, in turn, trapped 
by the light intensity profile of a laser tweezer. In this 
case the dielectric bead is free to rotate within the optical 
trap. On the other hand, attaching an end of a polymer 
to a super-paramagnetic bead, one can use magnetic field 
gradients to trap the polymer using a magnetic tweezer 
setup. In this case one can rotate the bead while holding 
it fixed in position by changing the direction of the exter- 
nal magnetic field. In the AFM experiments one end of 
a polymer is trapped by a functionalized tip of an AFM 
cantilever. The two distinct procedures which can be fol- 
lowed to measure force-extension are: (a) Both the ends 
of the polymer are held via the laser or magnetic tweezers 
or the AFMs. (b) One end of the polymer is attached to 
a substrate such that the position and orientation of this 
end is fixed while the other end is trapped via a laser or 
magnetic tweezer or an AFM cantilever. 

While the optical tweezers allow free rotation of di- 
electric beads within the trap, thereby, allowing free ori- 
entations of the polymer end, the magnetic tweezers fix 
the orientation of the ends and one can study the de- 
pendence of polymer properties on end-orientations by 
controlled change of the direction of external magnetic 
field. In this paper, we call this fixing of orientation of 
an end of a polymer as grafting. By changing the trap- 
ping potential from stiff to soft t rap one can go from the 
Helmholtz to the Gibbs ensemble[l3|. Before we proceed, 
let us first elaborate on how to fix the ensemble of a me- 
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chanical measurement [T^ In the simplest case we 
can assume that one end of the polymer is trapped in a 
harmonic well, V{z) — C{z — zo)^/2 with (0,0, zq) be- 
ing the position of the potential-minimum. The polymer 
end will undergo continuous thermal motion. One can 
use a feedback circuit to shift zq to force back the fluc- 
tuating polymer end to its original position. This will 
ensure a Helmholtz ensemble. This can also be achieved 
by taking C — > cx). On the other hand, one can use a 
feedback circuit to fix the force —C{z — zq) by varying 
C depending on the position z of the polymer end. This 
will ensure a Gibbs ensemble. This can also be achieved 
by taking a vanishingiy soft (C — *■ 0) trap to infinitely 
large distance {zq — » oo) such that within the length scale 
of fluctuation the polymer end feels a constant slope of 
the parabolic potential. Surely, in experiments, using a 
feedback circuit is easier to implement a particular en- 
semble. However, the other procedure is mathematically 
well defined and one can seek recourse of it to show that 
the partition function of two ensembles are related by a 
Laplace transform [l^ . This relation does not depend on 
the choice of the Hamiltonian for a polymer. An exact 
relation between the two ensembles for worm like chain 
(WLC) model is shown in SecHIl 

From the above discussion on possible experiments, it 
is clear that there can be three possibilities of boundary 
conditions in terms of orientations: (a) Free end: Both 
the ends of a polymer can remain free to rotate [l^ [l^ . 
(b) One end grafted: Orientation of one end is fixed and 
the other can take all possible orientations [l7| . (c) Both 
ends grafted: Orientations of both the ends are kept 
fixed. Thus, in experiments, one can have two possible 
ensembles and three possible boundary conditions. Wc 
restrict ourselves to the WLC polymers embedded in two 
dimensions (2D). We investigate the probability distribu- 
tion, free energy profile and force extension relation for 
each of these cases in this paper. We shall see that the 
properties of a semiflexible polymer depend both on the 
choice of the ensemble and the boundary condition. Note 
that, there can be other possibilities of boundary condi- 
tions e.g. orientation at one end of a polymer can be free 
to rotate on a half-sphere However, in this paper we 
focus on the three possible boundary orientations listed 
above. 

The WLC model is a simple coarse grained way to cap- 
ture bending rigidity of an unstretchablc polymer p^ . [20| 
embedded in a thermal environment. Recent siiigle 
molecule experiments in biological physics [1, |^, |^, |T3| 
renewed interest in this old model of polymer physics. It 
was successfully employed [2l|, [22j] to model data of force- 
extension experiments [1] on dsDNA. Mechanical proper- 
ties of giant muscle protein titin [H, [23| , polysaccharide 
dextrane [12, 24] and single molecule of xanthane [2^ 
were also explained using the WLC model. Due to the in- 
extensibility constraint, the WLC model is hard to tract 
analytically except for in the two limits of flexible chain 
{t — !■ oo) and rigid rod {t 0) , about which perturbative 
calculations have been done [la, [23, [29|] . A key quan- 



tity that describes statistical property of such polymers 
is the distribution of end-to-end separation. Numerical 
simulations to obtain radial distribution function for dif- 
ferent values of t have been reported along with a se- 
ries expansion valid in the small t limit [30,]. Mean- field 
treatments to incorporate the inextensibilty in an ap- 
proximate way have also been reported [Sll [3^ . In an 
earlier studv|16j we investigated the free energy profile 
of a semiflexible polymer whose ends were free to rotate 
in the constant extension ensemble and in the stiffness 
regime of 1 < t < 10. This work predicted that a clear 
qualitative signature of semiflcxibility would be a non- 
monotonic force extension for stiffnesses around t ~ 4 
in the Helmholtz ensemble. This comes from the multi- 
modality of probability distribution of end to end sepa- 
ration. However, this non-monotonicity is absent in the 
Gibbs ensemble [16|. Multiple maxima in the probabil- 
ity distribution of end to end separation was due to a 
competition between entropy, that prefers a maximum 
near zero separation, and energy, that likes an extended 
polymer. A series of later studies [l^, HI, [H, [11] used an- 
alytic techniques to understand the end to end distribu- 
tion at all stiffnesses including the stiffness regime where 
multimodality was observed. Recently, multimodality is 
found in transverse fluctuations of a grafted polymer us- 
ing simulations T7] and approximate theorv[36l [STj. A 
Greens function technique has been developed that takes 
into account the orientations of the polymer ends [sst . 
The impact of the specific boundary conditions and the 
comparable length scales of a dsDNA and the beads to 
which it is attached in typical force-extension measure- 
ments have been identified in another recent study [l^. 
The WLC model has also been extended to study statis- 
tics of end to end separation and loop formation prob- 
ability in dsDN A[3S| a. nd to incorporate twist degree of 
freedom [H [H [tfll . 

The construction of this paper is as follows. In SecHIl 
we present a theoretical technique for exact calculation 
of the WLC model via a mapping to a quantum parti- 
cle moving on the surface of an unit sphere. This tech- 
nique incorporates all the possible end orientations and 
predicts results in both the Helmholtz and the Gibbs en- 
sembles. In Sec lIIII we discuss the different discretized 
versions of the WLC model and the Monte-Carlo (MC) 
simulation procedures followed in this work. In Sec lIVI 
we present all the results of probability distributions and 
force-extensions etc. obtained from theory and simula- 
tions. Then, in Sec(Vl we summarize our results and 
conclude with some discussions. 



II. THEORY 

In the WLC model a polymer is taken as a continuous 
curve denoted by a d-dimensional vector r(s) where s is 
a distance measured over the contour of the curve from 
one of its ends. This curve has a bending rigidity and 



3 



thus the Hamiltonian is given by 



(1) 



where t{s) = dr{s)/ds is the tangent vector and the poly- 
mer is inextensible i.e. P = 1, /3 is the inverse tem- 
perature. Persistence length is a measure of the dis- 
tance up to which the consecutive tangent vectors on 
the contour do not bend appreciably and is defined by 
(£(s).t(0)) — exp(— s/A). The bending rigidity n is re- 
lated to persistence length A via n = [d— l)A/2. 

In this section we present a theoretical method to solve 
the WLC model to any desired accuracvfisl [4^ for both 
the Helmholtz and the Gibbs ensembles and all the three 
possible boundary orientations over the entire range of 
stiffness parameter t. We first present the method for a 
free polymer [l^. Then we extend it to calculate proper- 
ties of grafted [ one/both end(s) ] polymers. 

The partition function of a WLC polymer in the 
Helmholtz ensemble is Z{r) = ^^exp(— /?7i) where c 
denotes a sum over all possible configurations of the 
polymer that are consistent with the inextensibility con- 
straint. The probability distribution of the end to end 
vector becomes, P{f) = Z{r)/ dfZ{r) = MZ{r). If 
the tangent vectors of the two ends of a polymer are held 
fixed at ti and t/, the probability distribution of end to 
end vector in constant extension ensemble can be written 
in path integral notation as 



P(f) = U / V[i{s)] exp {-pn) X 5'^ 



where I?[i(s)] denotes integration over all possible paths 
in tangent vector space from the tangent at one end ti 
to the tangent at the other end tf. In c?-dimensions r = 
{fiT^^T ■ ■ , fd)- Recently a path integral Greens function 
formulation has been developed [l^ to evaluate the end 
to end distribution for a free polymer in 3D. We closely 
follow that method and generalize it to obtain results 
for various orientation constraints on polymer ends. In 
particular we focus on polymers living in a 2D embedding 
space. 

The integrated (projected) probability distribution is 
given by, 



P^{x) = / drP{r}S{ri-x). 



(3) 



We define the generating function of Px{x) via a Laplace 
transform, 



Pif) = / rfa:cxp(/a;/A)P,(x) 



(4) 



where / is the force in units of ksT/X i.e. / = FX/ksT 
applied along the a;-axis. Again, the partition function in 
the Gibbs ensemble, Zif) = dfexpif.f/X)Z{f)^. 



This immediately gives, = 1/Z{0). We show that 
Z{0) is a constant which depends on the constraints on 
end orientations. Eq|4] gives. 



P{f) = Af / V[i{s)]e 



U / V[i{T')]e 



(5) 

The last step is obtained by replacing r' = s/A and using 
the identities k = {d — l)A/2 and t = L/X. Note that, 
P{f), is the partition function, apart from a multiplica- 
tive constant, in the Gibbs ensemble where t behaves like 
an inverse temperature such that the Gibbs free energy 
can be written as G(/) = — 1/i lnP(/). Now considering 
t' as imaginary time and replacing t — —it' one gets, 
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with the identification of C = 

Lagrangian, P{f) [ — Z{f) /Z{Q)] in the above expression 
is the path integral representation for the propagator of 
a quantum particle, on the surface of a d- dimensional 
sphere, that takes a state to \if). In Schrodinger 
picture this can be written as the inner product of a state 
\ti) and another state evolved by imaginary time —it, 



ids (2) Zif) = {U\eM-^Hi-it))\if) = {U\exp{-tH)\if), (7) 



where H is the Hamiltonian operator corresponding to 
the Lagrangian C. 

Once P{f) = Z{f)/Z{0) is calculated, performing an 
inverse Laplace transform one can obtain the projected 
probability distribution P^ (x) . Eq|3] can be written as, 



(8) 



P{.f) = / dvxexp{tf Vx)pxivx) 



where Vx — x/L and Px{vx) = LPx{x) is a scaling re- 
lation. Note that the Helmholtz free energy is given by 
^x{vx) — —{^/t)\npx{vx)- Thus EqlUgives the relation 
between the Helmholtz and the Gibbs ensemble for finite 
chain (finite t). 



exp[-iG(/)] 



J dvxexp{tf Vx)exp[-tJ^x{vx)]- 



In thermodynamic limit of i ^ cx), a steepest descent ap- 
proximation of the above integral relation gives G{f) = 
^x{vx) — fvx, the well known Legendre transform rela- 
tion. Identifying —iu = tf one can define Fourier trans- 
form relations, Px{u) = J^j^px{vx) Bxp{—iuvx)dvx and 
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Pxivx) ^ -z- / dupxiu) exY){iuVx) 



(9) 
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FIG. 1: (Color online) For a semiflexible polymer in 2D hav- 
ing its ends free to rotate Px{vx) ( = Py{vy) ) is plotted 
at stiffness parameter t = 2. The points are collected from 
Monte-Carlo simulation in freely rotating chain model (see 
Sec lIIip . The line is calculated from theory (see Sec|ll]). The 
theory shows excellent agreement with simulation. It clearly 
shows bimodality via two maxima in integrated probability 
distribution at the two near complete extensions. 



such that P{f) = Px{u = ift) and the inverse Fourier 
transform can be written as an inverse Laplace transform, 



Px(Vx) t — 
2m 



dfP{f) eM~tfvx) 



(10) 



The simplest way to obtain px{vx), numerically, is to 
replace / — —iu/t in the expression for P{f) to obtain 
Px(u) and evaluate the inverse Fourier transform (Eql9]). 

Up to this point everything has been treated in d- em- 
bedding dimensions. Experiments on single polymer can 
be performed in three dimensions (3D) as well as in two 
dimensions (2D). In 3D, polymers are left inside a solu- 
tion whereas one can float the polymer on a liquid film 
to measure its properties in 2D 0]. However, polymers 
embedded in 2D are more interesting because of the fol- 
lowing reason. In a free polymer whose end orientations 
are free to rotate, the system is spherically symmetric 
and thus the probability distribution of end to end vec- 
tor P{f) = P{r) where r = |r|. For this system it was 
shown that in the Helmholtz ensemble in 3D|15|. 



, I dpx 
V{vx) = -r- 

2'KVx dVx 



(11) 



where p{v = r/L) = L'^P{r) is the probability distribu- 
tion of end to end distance scaled by the contour length 
L. In presence of the spherical symmetry of a free WLC 
polymer, this distribution gives the Helmholtz free en- 
ergy T{v) — — (1/t) lnp(-y)[lf^]. P{r) is related to the 
radial distribution function S{r) via S{r) — Cdr'^^^P{r) 
where Cd is the area of a d- dimensional unit sphere. 
Since p{v) is a probability distribution, p{vx) > and 
therefore dpx/dvx < for Vx > thus ruling out mul- 
tiple peaks in Px{vx) and showing that Pxivx) will 



have a single maximum at w^. = for all values of stiff- 
ness parameter t. No such simple relation exists between 
p{vx) and Px(vx) in 2D. The two dimensional WLC poly- 
mer having its ends free to rotate may show more than 
one maximum in Pxi^x) and therefore non-monotonicity 
in force-extension. Indeed our calculation and simulation 
(see Sec lIII[) does show multiple maxima in projected dis- 
tribution Pxivx) (Fig[T]). This is a curious difference be- 
tween semiflexible polymers in 2D and 3D. Because of 
this and the fact that experiments in 2D are possible 
in this work we focus on the 2D WLC polymers. 

We have already given a general form of Z{f ) (Eq[7|) 
which depends on the dimensionality d of the embed- 
ding space. For d = 2, one can assume t = (cos 6*, sin 0), 
leading to £ = {1/4 -I- fcos6}. This automatically 
maintains the inextensibility constraint P — 1. The an- 
gular momentum pg = ^ ^ 9/2 and thus the corre- 
sponding Hamiltonian H — Opg — C = Pg — fcosO. In 
planar polar coordinates, replacing pg — s- —i-§g one ob- 
tains the corresponding quantum Hamiltonian operator, 
^ = ~"§ffi ~ fcosO. In this representation of tangent 
vectors, 



z{f) = {e.,\cM-tH)\ef) 

= Y.'t'nmK'{0f){n\eM-tH)\n') 



(12) 



where 0„(6') = {n\6). If external force is applied along a;- 
direction as in EqHl H = Hq + Hj = --§^- fcosO. Thus 
the total Hamiltonian H denotes a rigid rotor [Hq = 
— -§gi) in presence of a constant external field {Hj = 
—fcos6). The eigenvalues of Hq are En = and the 
complete set of orthonormal eigenfunctions are given by 
(j^niO) = exp(m0)/\/27r where n = 0, ±1, ±2, . . . , ±oo. In 
this basis {n\Hi\n') = -(//2)(5„/,„+i + 5n\n-i)- There- 
fore, {n\H\n') = n25„/_„- (//2)((5„',„+i + (5„',„-i). If the 
external force were applied in y- direction Hi = — f smO 
and {n\H\n') = n^(5„',„ - (//2i)((5„',„+i - Sn'.n-i)- 
(n| exp(— i77)|n') can be calculated by exponentiating the 
matrix {n\H\n') . Thus one can find Z{f) and hence P{f) 
and px{vx)- 

Note that the above formalism can be easily extended 
to find the end to end vector probability distribution 
p{vx,Vy). A Laplace transform of P{r) is P{f) — 
dfeKp(r.f/X)P(r). In a similar manner as above one 
can show that P{f) = Z{f)/Z{Q) with Z{f) given by 

EqlT2]with H — — ^ — fx cos 9 — fy sinO. Thus, using 
an inverse Laplace transform one can find P(f) and hence 

P{Vx,Vy). 

A. Free polymer 

For a polymer which has both its ends free to ro- 
tate, integrating EqlT^] over all possible initial and fi- 
nal tangent vectors in rigid rotor basis one gets, Z{f) = 
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2tt (0| exp(-<ff)|0), Z(0) = 2tt and hence 



III. SIMULATION 



Hf) 



\eM-tHm 



(13) 



This means that P{f ) is given by the (0, 0)-th element of 
the matrix (n| exp(— Thus, if the external force 
/ is applied in x- direction, remembering Px{u) — P{f ~ 
~iu/t) one can calculate the inverse Fourier transform 
(EqlH]) to obtain Px{vx)- In this case, due to spherical 
symmetry of a polymer whose ends are free to rotate, 

Px{Vx) = PyiVy)- 



B. One end grafted 

This symmetry breaks down immediately if one end of 
the polymer is fixed to a specific direction, namely along 
the X-axis i.e. 6i — 0. Then in Eq. ()12p integrating over 
all possible Of and leaving 9i = one obtains Z{f) = 
^„(n| exp(— t_ff)|0) in the rigid- rotor basis. Note, for 
this case Z{0) ~ 1 and therefore 

Pif) = J2{n\cM-tH)\0) . (14) 



C. Both ends grafted 

Two ends of a polymer can be grafted in infinitely 
different ways. Let us fix the orientation of one 
end along x- direction {9i = 0) and the other end 
along any direction 9f. Then Eq. (fT2|) gives 2nZ{f) — 

E„,„'e™''^(n|exp(-tif)|n'), 2^Z(0) = ^„e™«/-*"' 
and hence 



n ^ 



(15) 



If the external force is in x- direction, the Laplace 
transform of Z{f), defined in the way described above, 
gives the projected probability distribution in x- direc- 
tion, Px{vx)- On the other hand, if the external force is 
in y- direction, the Laplace transform of Z(f) gives the 
projected probability distribution in y- direction Pyivy), 
the distribution of transverse fluctuation while one end 
of the polymer is grafted in x- direction. 

All the relations derived so far are exact. Since the 
calculation of an infinite order matrix (n| exp(— t_ff)|n') 
is not feasible, we calculate it numerically [45] by truncat- 
ing up to an order iVd, that controls the accuracy, limited 
only by computational power. Unless otherwise stated, 
we use Nd — 11 which already gives very good agreement 
with simulated data (see Fig[l]and Sec lIVp . The inverse 
Laplace transforms to obtain end to end probability dis- 
tributions from P{f)s are also done numerically. 



In this section, we introduce two discretized models 
that we use to simulate semiflexible polymers. Both of 
these are derived from the WLC model which has been 
used for our theoretical treatment in SecHIl After in- 
troducing the discretized models we show how to impose 
the various boundary conditions on end orientations. We 
perform Monte-Carlo (MC) simulations of these models 
to obtain probability distributions in the Helmholtz en- 
semble. 

One discretized version of the Fokker-Plank equation 
corresponding to the WLC model is the freely rotating 
chain (FRC) model [H, [13 • In the FRC model, one con- 
siders a polymer as a random walk of N steps each of 
length b = L/N with one step memory, such that, suc- 
cessive steps are constrained to be at an fixed angle 9 
with A = 2h/9^ . The continuum WLC model is obtained 
in the limit 0, & ^ 0, TV — *■ oo keeping A and L finite. To 
simulate a polymer with ends free to rotate a large num- 
ber of configurations are generated with first step taken 
in any random direction. Whereas if one chooses the first 
step to be in some specific direction, this will simulate a 
polymer with one end grafted in that direction. 

A straight forward discretization of the Hamiltonian 
in Eq[T] in 3D (2D) is an Id Heisenberg (classical XY) 
model: 



f3n 



N 



N 



with a nearest neighbor coupling J = n/h between 'spins' 
ti. We have ignored a constant term in energy. The ap- 
propriate continuum limit is recovered for & — > 0, J ^ oo 
with Jb = K finite. In this model grafting is simulated by 
fixing end spins on the ID chain. If an end is free then 
the end spin takes up any orientation that are allowed 
by the energy and entropy. In this model, by fixing the 
two end-spins, one can easily simulate a polymer with 
both its ends grafted in some fixed orientations. We fol- 
low the normal Metropolis algorithm ^4^ to perform MC 
simulation in this model. 

We restrict ourselves to two dimensions. In the FRC 
model simulations we have used a chain length of = 
10'^ and generated around 10^ configurations. This sim- 
ulation does not require equilibration run. Therefore all 
the 10^ configurations were used for data collection. In 
the XY model we have simulated A^ = 50 spins and equi- 
librated over 10^ MC steps. A further 10^ configurations 
were generated to collect data. We have averaged over 
10'^ initial configurations, each of which were randomly 
chosen from nearly minimum energy configurations that 
conform with the boundary conditions. Increasing A^, in 
both the models of simulation, do not change the aver- 
aged data. As a check on the numerics, we compared 
simulation evaluation of (r^) and (r**) with their exact 
results[l^ |3^ to obtain agreement within around 0.5%. 

Notice that in simulating the FRC model one performs 
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FIG. 2: (Color online) None of the force-extension curves, ob- 
tained in the Gibbs ensemble, including that at f = 3.33 show 
non- monotonic behavior unlike in the Helmholtz ensemble 
[l^ . Forces are expressed in units of fcsT/A i.e. / — FX/kBT. 

random walk with fixed angle between consecutive steps 
and does not require to equilibriate. Thus one uses all 
the simulated configurations for data collection. On the 
other hand, in simulating XY-model one has to perform 
equilibration runs over a large number of steps and an 
averaging over many initial configurations is required. 
Another important difference between the two simulation 
methods is that, in the XY-model simulation, in each MC 
step one has to calculate a time consuming exponential 
of change in energy, whereas, no such exponential calcu- 
lation is required in simulating the FRC model. Thus 
simulating the FRC model is clearly much faster, com- 
putationally. However, implementing the fixed boundary 
orientations at both the ends of a polymer is much easier 
in the XY-model. 



IV. RESULTS 

Once all these theoretical and simulation tools are 
available, we apply them to bring out the statistical and 
mechanical properties of a semiflexible polymer. We have 
three different boundary conditions depending on the 
orientational constraints on the polymer ends and two 
different ensembles. For each case we look at the vari- 
ous probability densities, ensemble dependence of force- 
extension etc. For the case of a polymer with both ends 
grafted we find that the properties depend on the relative 
orientation of the two ends. 



A. Free polymer 

The Helmholtz ensemble: We employ the theory as de- 
scribed in Sec|TT]to calculate Px{vx) and Py{vy) for a poly- 
mer with both its ends free to rotate. We compare the 



probability distributions obtained at stiffness parameter 
t = 2 with that obtained from MC simulation fSec lIIip 
using the FRC model (see FiglT]). This shows excellent 
agreement between theory and simulation. For a free 
polymer Px{vx) and Py{vy) are same due to the spheri- 
cal symmetry. Note that J-{vx) = — (1/t) hipx{vx) would 
give a non- monotonic force-extension {fx)-Vx due to the 
multimodality in Px{vx) (FigH]) via {fx) = {dT/dvx). 
The force-extension obtained from the projected proba- 
bility distribution (ua;) corresponds to the experimental 
scenario in which the external potential traps the poly- 
mer end only in the x- direction and the polymer-end 
is free in y. In general, if the external potential traps 
the polymer-end in dr dimensions {dr < d) then a dr di- 
mensional projection {[d — d^] dimensional integration ) 
of the probability distribution of end to end vector p{v) 
gives the appropriate free energy and decides the force- 
extension relation. On the other hand, if the trapping 
potential holds a polymer-end in all the d-dimensions, as 
is usually done in most force-extension experiments, only 
the end to end vector distribution p{v) gives the appro- 
priate Helmholtz free energy that can predict the force- 
extension behavior in the Helmholtz ensemble. This un- 
derstanding is general and does not depend on the spe- 
cific orientational boundary conditions or the dimension- 
ality d of embedding space. This is important to keep in 
mind while analyzing experimental data. In experiments 
that use the laser tweezers to trap polymer ends in d- 
dimensions, ends remain free to rotate and the relevant 
Helmholtz free energy is obtained from p{v). Ref.flGj 
predicted multiple minima in this free energy leading to 
non-monotonic force-extension in such experiments. 

The Gibbs ensemble: We have already mentioned that 
the non- monotonic nature of force-extension, a strong 
qualitative signature of semiflexibility, is observable only 
in the Helmholtz ensemble and not in the Gibbs en- 
semble [l^. In the Gibbs ensemble, the averaged ex- 
tension comes out to be (v) = ~{dG/df ) and the re- 
sponse d{v)/df — t[{v'^) — {v)"^] > 0. Similar relation 
for response function does not exist in the Helmholtz 
ensemble. Therefore, the force-extension in the Gibbs 
ensemble has to be monotonic (Fig[2]) in contrast to the 
Helmholtz ensemble. For a polymer with its ends free 
to rotate, the force extension relations, that have been 
calculated from theory, at various t are shown in Figl^l 
For small forces the polymer shows linear response. At 
large and positive force polymer goes to fully extended 
limit beyond which, the inextensibility constraint pre- 
vents further extention. It is possible to do perturbative 
analysis of P{f) — (0| exp(— ii7)|0) in the two extreme 
limits of small and larg e forces to obtain the asymp- 
totic force-extensions |47j . In the small force limit, f cos9 
may be treated as a perturbation about the rigid ro- 
tor hamiltonian Hg = —d^/d9^. Thus keeping upto 
the second order correction to eigen-values we obtain 
Eq = — /^/2. Within this perturbative approximation 
P{f ) = exp{-tEo) and therefore G(/) -1/t lnP(/) = 
— /^/2. Thus the force extension relation in this limit 
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FIG. 3: (Color online) The simulation data for p^iv^) and 
Py{vy) from the FRC model and the XY model are compared 
with their theoretical estimates. Simulations and calculations 
were done at t — 2 for a polymer with one end grafted in x- 
direction. 



is (v) — "dG/df — /. On the other hand, for large 
forces one can expand the term cos 61 ~ 1 — 9^/2 and 
write H = -tf + tHo where Ho = -d^/d9^ + (1/2) f 9^ 
is the harmonic oscillator Hamiltonian. In the harmonic 
oscillator basis, the ground state eigenvalue Eq = \/,f /2 
and thus the ground state energy corresponding to H is 
^ / + \fTJ^- Therefore, in a similar manner as in above, 

P(/) =cxp(-<V772 + t/) andG(/) = ^[Jl2- ! . Thus, 
for large forces, the force-extension relation comes out to 
be (v) = —dG/df = 1 — 1/(2^/2/) which can be inverted 
to get the relation, / — 1/[8(1 — (w))^]- AH the curves 
f({v)) in Figl2] falls on to / = (v) at / — > limit and to 
1/[8{1 - in the / ^ oo limit. 



B. Grafted polymer: One end 

The Helmholtz ensemble: Let us compare our theoreti- 
cal and simulation estimate oipx{vx) sjhA py(vy) dXt — 2 
(Figl3]) for a semiflexible polymer with one end grafted 
in X- direction. The excellent agreement validates both 
our theory and the simulation techniques. In Px{vx), the 
peak in near complete extension along positive x is due 
to the coupling of the end orientation towards this di- 
rection with large bending energy (also see Fig fT^ . We 
then explore, in detail, the transverse fluctuation Py{vy) 
of this system for different t (FigH]). At large t{= 10), 
Py{vy) has single maximum at Vy = 0. At such low stiff- 
nesses entropy takes over energy contributions. Number 
of possible configurations and thus entropy gains if end 
to end separation remains close to zero. This gives rise to 
the single central maximum. The emergence of multiple 
maxima at nonzero Vy, the multimodality, at larger stiff- 
ness {t = 2.8) is due to the entropy- energy competition. 
The central peak is due to the entropy driven Gaussian 
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FIG. 4: (Color online) For a polymer with one end grafted in 
X- direction, the integrated probability distribution Py{vy) is 
plotted at various stiffnesses t. At t = 4 there is a single max- 
imum at Vy = 0. Decreasing t we see at t = 2.8 emergence of 
two more peaks at nonzero Vy apart from the one at = 
(See inset). At t — 2 the central peak vanishes, the trimodal 
distribution becomes bimodal. The circles labeled LMF are 
data taken from Ref. at t — 2 and show excellent agree- 
ment with our theory. At t — 0.75 we see re-emergence of the 
central peak and trimodality in Py{vy) (See inset, Os are from 
our MC simulation in the FRC model at t — 0.75.) The lines 
are calculated from theory. 




FIG. 5: (Color online) The left panel shows the Helmholtz 
free energies J-{vx) and J-{vy) of a polymer at t — 2 and one 
end grafted in x- direction. The right panel shows the corre- 
sponding force-extensions in the Helmholtz ensemble. Both 
(fx)- Vx and (/a)- Vy show non- monotonicity and regions of 
negative slope. Free energies are expressed in units ot ksT 
and forces are expressed in units of ksT/X. 



behavior. The other two peaks emerge as entropy tries to 
fold the polymer and energy restricts the amount of bend- 
ing. Since bending in positive and negative y- directions 
are equally likely, the transverse fluctuation shows two 
new maxima near Vy = ±0.5 symmetrically positioned 
around Vy — 0. With further increase in stiffness {t — 2), 
the central entropic peak vanishes (also see Fig fT2)) and 
Py{vy) becomes bimodal with two maxima (Fig[3]). At 
even higher stiffness {t = 0.75) the central peak reap- 
pears, due to a higher bending energy. At t = 0.5 the 
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FIG. 6: (Color online) Average displacements along x- di- 
rection (vx) and y- direction (vy) as a function of transverse 
force (transverse to grafting direction a;) in constant force en- 
semble. Lines denote our theoretical calculation while points 
denote the MC simulation data taken from Ref. [l3]. Forces 
(Fy) are expressed in units of ksT/X, i.e. fy — FyX/ksT. 

distribution again becomes single peaked at Wy = as 
bending energy takes over entropy and the polymer be- 
comes more like a rigid rod. However, even at very high 
stiffness like t = 0.5 the single peaked distribution Py{vy) 
is quite broad underlining the influence of entropic fluctu- 
ations. Notice that we have plotted MC data taken from 
Ref.tlTi] for the XY model simulation at i = 2 (FigH]). 
This shows very good agreement with our theory. Infact 
all the simulated data from Ref. at different t show 
excellent agreement with our theoretical predictions. In 
the inset of FigHl we have magnified the multimodality 
at t = 2.8 and t = 0.75. We have also plotted our PRC 
model simulation data ed t — 0.75 and obtained very 
good agreement. 

At this point, it is instructive to look at the force 
extension behavior in the Helmholtz ensemble, the en- 
semble in which Py{vy) and Px{vx) have been calcu- 
lated above. In it the extension Vx [vy\ is held con- 



stant and the corresponding average force in x- [y-] di- 
rection is found from the relation {fx) = dP{vx) / dvx 
( or {fy) = d!F{vy)/dvy ). Notice that, when Vx [vy] is 
held constant, Vy [vx\ remains free. This can be achieved 
using a trapping potential constant in Vy [vx\ and trap- 
ping the polymer end in Vx [vy\. In FigO we show the 
Helmholtz free energies J-(vx) = — (l/t)ln Px{vx) and 
^{vy) — — (l/t)ln Py(vy) and the corresponding force 
extension curves in constant extension ensemble. Note 
that unlike the monotonicity obtained in {Vy)-fy curve 
(FiglB]) in the Gibbs ensemble, the {fy)'Vy curve in FigO 
clearly shows non-monotonicity, a signature of semiflexi- 
bility in the Helmholtz ensemble. 

The Gibbs ensemble: From our theory we can also ex- 
plore the transverse response of a polymer which has one 
of its ends grafted and a constant force is applied to the 
other end in a direction transverse to the grafting direc- 
tion. Assume that the grafting direction is x and a force 
fy is applied in y- direction to study the transverse re- 
sponse. A linear response theory was proposed earlier (soj 
to tackle this question. Our theory can predict the effect 
of externally applied force fy of arbitrary magnitude on 
the average positions {vx) and {vy). As the force is ap- 
plied in y-direction i.e. / = yfy, we have Hj = — fy sinO. 
Because one end of the polymer is grafted in x- direc- 
tion we use {n\Hi\n') = -{fy/2i){dn',n+i - Sn',n~i) to 
evaluate Z{fy), whereas to calculate {vx) = —{dG/dfx) 
[ or, {vy) ~ —{dG/dfy) ], we introduce a small perturb- 
ing force dfx [ or, Sfy ] in the Hamiltonian matrix to 
obtain the partial derivatives. Thus we obtain the corre- 
sponding force-extensions shown in FiglSj As the grafted 
end is oriented in x- direction, we expect, in absence of 
any external force, {vx) will be maximum and will keep 
on reducing due to the bending of the other end gen- 
erated by the external force fy imposed in y- direction. 
Thus {vx) is expected to be independent of the sign of 
fy. Similarly, {vy) should follow the direction of exter- 
nal force and therefore is expected to carry the same 
sign as fy. Fig [H] verifies these expectations and shows 
very good agreement between our theory and simulated 
data taken from Ref. [l7j . It is interesting to note that, in 
the Helmholtz ensemble, the multimodality in probability 
distribution predicts non-monotonicity in force-extension 
relation. However, as expected, this non-monotonicity 
does not survive in the Gibbs ensemble. 



C. Grafted polymer: Both ends 

The Helmholtz ensemble: Let us first fix the orienta- 
tions of the polymer at both its ends along x- axis and 
compare Pxi^x) and Py{vy) obtained from our XY model 
simulation and our theory (Fig[7]). The very good agree- 
ment validates both our theory and simulation. Then, 
we go on to explore the properties of this system using 
the theory developed in Sec|TT]-C. Let us fix the orienta- 
tion at one end in x- direction (9i = 0) and that in the 
other end (Of) can be varied to study the change in trans- 




FIG. 7: (Color online) The simulation data for px{vx) and 
Py{vy) from the XY model simulations of a WLC polymer are 
compared with their theoretical estimates. Simulations and 
calculations were done ai t — 2 for a WLC polymer with both 
its ends grafted in x- direction. 



verse fluctuation Py{vy). To begin with, let us find Py{vy) 
for different stiffness parameters t with 6f = (Fig[5]). 
The height of the central peak shows non-monotonicity 
- with increase in t from t — I the height of the central 
peak first decreases up to t = 2 and then eventually it 
increases again. The initial decrease in peak height is 
due to the fact that with increase in t, i.e. with lower- 
ing in stiffness, the other end of the polymer (relative to 
the first end) starts to sweep larger distances from the 
X- axis. With further increase in t (t = 4), the height of 
the maximum increases (also see Fig fT^ . From FigfT^ 
notice that at t = 4 multimodality appears in the distri- 
bution of end to end vector. The new entropic maximum 
at Tj = contributes towards increasing the peak height 
in Py{vy) at Vy — 0. Though, in p{vx,Vy) multimodality 
is present (see Fig fT^ a.t t — 4, after integration over 
probability weights along x- direction the projected dis- 
tribution Py{vy) becomes unimodal. Thus multimodality 
in the probability distribution of end to end vector does 
not guarantee multimodality in projected probability dis- 
tributions. 

To see the impact of change in relative angle of graft- 
ing, we fix one end along x- axis and rotate the orienta- 
tion of the other end and find out the transverse fluctu- 
ation ^,^(1;,^) at t = 4 (FigED. At 61 = 6*/ the fluctu- 
ation is unimodal with the maximum at Uj, = 0. With 
increase in 9 the orientation of the other end rotates from 
positive x- axis towards positive y- axis. Energetically 
the polymer gains the most, if it bends along the perime- 
ter of a circle. Therefore, energetically, at any 9, the peak 
oipy{vy) would like to be at Vy ~ { 1 — cos6' )/9. Thus at 
9 = 0, 7r/4, 7r/2, 37r/4, 77r/8, TT the peak oi py{vy) should 
be at Vy = 0, 0.37, 0.64, 0.72, 0.69, 0.64 respectively. FiglS] 
shows that the peak positions almost follow these values 
up to = 7r/2, above which entropic contributions dom- 




-1 -0.5 0.5 1 




V 



FIG. 8: (Color online) The upper panel shows Py{vy) for a 
polymer with both ends grafted along x- direction at various 
stiffness parameters t. They always show single maximum. 
In lower panel, Py{vy) is plotted for various relative angles 9 
between the orientations of the two ends ai t = 4. The inset 
magnifies the emergence of bimodality at ^ = 7r/4. 



inate to bring down the peak positions to lesser Vy with 
respect to that attained at 6* = 7r/2. However, entropy 
always play a crucial role, e.g. at 6* = 7r/4, Py(vy) shows a 
double peak around Vy — 0.37. At = tt the two ends of 
the polymer are kept anti- parallel. Notice that, as = tt 
and 9 = — tt are physically same, aX 9 — tt, energetically, 
Vy — ±0.64 are equally likely. Entropy would like the 
two ends to bend to Vy = 0. Competition between en- 
ergy and entropy leads to almost a constant distribution 
up to \vy\ ^ 0.5. The behavior of Py{vy) for — tt < 9 < 
is mirror symmetric about Vy = with respect to the 
behavior of Py{vy) in the region < 9 < n. 

The Gibbs ensemble: We then work in the constant 
force ensemble by applying a force / ~ yfy on an end 
oriented along any direction 9 to x-axis while the other 
end is oriented along x- direction. We find out the cor- 
responding responses, {vx)-fy and {vy)-fy to this force 
(FigjS]) in the similar manner as has been done in the 
last subsection for the case of a polymer with one end 
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FIG. 9: (Color online) Average displacements {vx) (upper 
panel) and {vy) (lower panel) as a function of a force fy for 
a polymer having one end grafted along x- direction and the 
other in an angle 6 to the x- direction. Forces are expressed in 
units of kBT/X. All the force-extension curves are obtained 
at t = 1. 



FIG. 10: (Color online) The upper panel shows the variation 
of (vx) as a function of 8 and the lower panel shows the varia- 
tion of (vy) as a function of 6. (vx) and (vy) are calculated for 
stiffness parameters t — 1,2,3,4,5, 10. The thick solid line, 
in both the plots, show the expected behavior coming from 
energetics ignoring the entropy. 



grafted. If 6* = 0, the force extensions carry the same 
quaUtative features as for a single end grafted polymer 
at all t (see 9 — curves for t = 1 in FigE]). Therefore, 
instead of showing the t dependence of force-extension 
behavior, we show the 6 dependence of force extensions 
a,t t — 1. The peak in {vx)-fy curve shifts to fy < as 
is increased up to 7r/2 above which it again shifts back 
towards fy — 0. With increase in 6, {vx) decreases, as 
with these boundary orientations the polymer is forced 
to close in x- and open up in y- direction. However, for 
9 — > ±7r entropy likes (vy) 0. For 6 < tt/2 small 
negative fy leads to unfolding thereby increasing (vx). 
Whereas for 6 > tt/2 the effect of negative force is oppo- 
site - it helps the polymer to get folded to reduce (vx). 
At = TT, (vx) always remains zero. The responses for 
negative 9 are reflection symmetric about fy = 0. The 
folding behavior is also apparent from {vy)-fy curves. Up 
to 61 = 7r/2 the response shifts towards positive (vy) as the 
polymer likes to open up in y- direction due to the bend- 
ing energy cost. However, for large 9 entropy wins and 
at = TT, {vy)-fy curve, again, goes through origin. The 



elastic constant dfy/d{vy) near fy — (linear response) 
is larger at 6* = as compared to at 6* = tt; i.e. the trans- 
verse response of a semiflexible polymer with parallel end 
orientations is more rigid than with anti-parallel end ori- 
entations. To see the impact of the change in relative 
angle 9, in detail, we calculate (vx) and (vy) as we vary 9 
(FigllOp keeping external force at zero. Bending energy 
would like {vx) = sm9/9 and {vy) = (1 — cos9)/9. Note 
that at 6* ^ 0, energetically, {vx) — > 1 and (vy) — > 9/2. 
Again, at 9 ±7r bending energy requires (vx) — > and 
(vy) ±2/7r though entropy likes (v) 0. Thus at 
small t, the approach of {vx)-9 curve to sn\9/9 is much 
better than approach of {vy)-9 to (1 — cos9)/9 (Fig fTU)) . 
It should be noted that the angle 9 in this study de- 
notes a relative angle of bending between the two end 
orientations of a WLC polymer. This should not be con- 
fused with the twist angle as in Ref. [41.]. In an earlier 
studyjlll the impact of changing 9 on the averaged root 
mean squared end to end vector has been obtained. In 
this section we have shown the impact of changing 9 on 
projected probability distribution, averaged end to end 
distance {{vx), (vy)) and force-extension relations. 
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FIG. 11: (Color online) The probability distribution of end 
to end distance 2-Kp{v) at stiffness t = 4 is plotted for the 
three different boundary conditions - (a) both ends free, (b) 
one end oriented in x- direction and the other kept free, (c) 
both ends oriented in x- direction. Radial distribution of first 
two cases are equal, whereas for the third case it is different. 
However, all the three curves show double maxima. 



D. Distribution of end to end vector 
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We now employ MC simulations to study some other 
aspects of probability distribution. We first examine 
the probability distribution of end to end distance p{v). 
It is clear from FiglTT] that grafting one end does not 
change the double maxima feature in p{v) at intermedi- 
ate values of stiffnesses (4 < t < 2). This is because 
the two cases are symmetry related; fixing orientation 
at one end only shifts the probability weight distributed 
over all possible angles at a given radial distance v to- 
wards the direction of the orientation. Though grafting 
both the ends change the distribution of end to end dis- 
tance, the double maxima feature persists and becomes 
more pronounced. We note that once one end of a poly- 
mer is grafted immediately the system loses its spherical 
symmetry, more so, since we restrict ourselves to semi- 
flexible regime. For a free polymer p(v) plays the role 
of a probability distribution of end to end vector and 
thus gives the Helmholtz free energy and force-extension 
behavior. Once the spherical symmetry is broken p{v) 
merely plays the role of a radial distribution function in 
terms of 2Trvp{v) and no longer remains relevant in pre- 
dicting the force-extension behavior. We have already 
seen that the projected probability distributions Pxivx) 
and py(vy) are very different for grafted polymers, though 
they are the same for free polymers that preserve spher- 
ical symmetry. 

The full statistics of the WLC polymers are encoded 
in the end to end vector distribution function piv^^Vy). 
To see the complete structure, we next obtain p{vx,Vy) 
from MC simulations in the PRC (for a free polymer or 
a polymer with one end grafted) and the XY model (for 



FIG. 12: (Color online) Density plot of p{vx,Vy). Color code: 
red (light) - high density, blue (dark) - low density. Left panels 
are for free polymers, middle panels are for polymers having 
one end grafted in x- direction and the right panels are for 
polymers having both ends grafted in x- direction. From top 
to bottom four panels denote increasing stiffness parameters 
t = 0.5,2,4,10 (decreasing stiffness). Note that the double 
maxima feature in p{vx,Vy) (one maximum near the centre 
and another near the rim) at t = 4 persists for all the three 
boundary conditions. 



a polymer with both ends grafted) and present them as 
two dimensional density plots. We compare p{vx,Vy) of 
a free polymer, a polymer with one end grafted and a 
polymer with both ends grafted (Fig fT^ . For definite- 
ness, we chose all the graftings, fixing of end orienta- 
tions, to be in the x- direction. We plot p{vx,Vy) over 
a range of stiffnesses {t = 0.5, 2, 4, 10). The distribution 
has finite values for v < 1 and is zero for t; > 1. This is 
due to the inextensibility constraint in the WLC model. 
In these density plots high probability is shown in red 
(light) and low in blue (dark) (Figdlj). At small stiff- 
ness [t = 10) p{vx,Vy) shows a single entropic peak at 
= for free polymer (Fig fT^ . This is slightly shifted 
towards the direction of end- orientations in grafted poly- 
mers. This shifted entropic peak slowly moves towards 
t/ = in the t — > oo limit. With increase in stiffness 
(t = 4), a new energy dominated probability peak ap- 
pears near the full extension limit, u = 1, of the polymer 
(Figll2p. This peak forms a circular ring for free poly- 



mers. For a grafted polymer, this new peak is aligned 
in the direction of grafting. The probability distribution 
p{vx,Vy) at i = 4 clearly shows two regions of proba- 
bility maxima, one near the zero extension and another 
near the full extension, for polymers with all kinds of 
boundary orientations - the free polymer, the polymer 
with one end grafted and the polymer with both ends 
grafted. Grafting of polymer ends, in a sense, enhances 
the effective stiffness. Therefore with increase in polymer 
stiffness (decrease in t) the multimodality sets in first in 
the polymer with both the ends grafted in the same di- 
rection near t = Q. At this t value the free polymer and 
the polymer with one end grafted show only the entropic 
peaks near v — Q. Near t = 5, the polymer with one 
end grafted starts to show multimodality. For free poly- 
mers, multimodality sets in only at an even higher stiff- 
ness close to t = 4. These behaviors are also borne out 
by the theory. The multimodality (two maxima) in prob- 
ability distribution of end to end vector seen for a free 
polymer at t = 4 (FigfT^ gives rise to the triple minima 
in free energy found in Ref.[16l]. In an earlier work[38l| 
it was shown that the crossover from flexible chain to 
rigid rod via multimodality in p robability distribution as 
obtained for a free polymer |16| persists even after graft- 
ing one end of the polymer. Here we have shown that 
this behavior persists even after grafting both the ends 
of a polymer. Certainly the detailed features of probabil- 
ity distribution of end to end vector would change with 
changing the relative angle of grafting at the two ends. 
At an even larger stiffness (t = 2), the entropic maximum 
near the centre {v = 0) disappears fFig fT^ . For the free 
polymer, one energy dominated maximum gets equally 
distributed over all angles. This way the system uses its 
spherical symmetry to gain in entropy. For grafted poly- 
mers, probability maximum near the full extension fans 
a finite solid angle around the direction of grafting. The 
distribution around the grafting direction is narrower for 
the polymer with both its ends grafted along the same 
direction. This is due to a larger coupling between graft- 
ing and bending stiffness. This fact is more pronounced 
inp(u^,wy) at t = 0.5 (FigHg). 

As mentioned earlier, in the Hclmholtz ensemble the 
free energy is given by T{vx.,Vy) = -(1/i) \Ti[p{vx,Vy)]. 
This free energy will give the force-extension behavior if 
the ends are trapped in 2D plane at some points (0, 0) and 
{vx, Vy). In Fig[T3]we plot this free energy profile JF(0, Vy) 
at <: = 4 and compare the three different boundary condi- 
tions. This plot clearly shows that triple minima in free 
energy [16] prevails even after grafting one or both ends of 
a semiflexible polymer. In terms of force-extension what 
this triple minima means? If we start off with end to end 
vector at (0, 0) and increase \vy\^ for small extensions the 
ends would experience an attractive force between them- 
selves. Beyond a limit {\vy \ ^ 0.5) the ends would repel 
each other to take the system to the other minima at non- 
zero Vy. At very large extension, again, they would expe- 
rience an attractive force, governed by the inextensibil- 
ity constraint. Thus, in force-extension experiments on a 
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FIG. 13: (Color online) At t = 4 free energy profile T{0,Vy) 
corresponding to the probability distributions shown in Fig |12l 
are plotted. This clearly shows that the triple minima feature 
in free energy for a polymer with both ends free persists even 
after grafting one or both ends of the polymer. 




FIG. 14: (Color online) Density plot of p{vx,Vy) obtained 
from theoretical calculations. Color code: red (light) means 
high probability and blue (dark) means low. All the plots are 
obtained for stiffness t = 4; (a) ends free to rotate, (b) one 
end grafted along x- direction, (c) both ends grafted along x- 
direction. 



polymer in constant extension ensemble, this multistabil- 
ity (non-monotonic force extension) at intermediate stiff- 
ness values should be measurable for all kind of boundary 
conditions. However the measurement would require av- 
eragin g oy er a large number of observations as indicated 
in Ref . [16| . At this point, it is interesting to notice that, 
for a polymer with one end grafted along x-direction the 
force-extension obtained from the slope of J^(vx = 0, Vy)- 
Vy curve gives the transverse response in constant exten- 
sion ensemble when the other end is constrained to be 
at a fixed v = {0,Vy). The behavior of the transverse 
response, evidently, would then also depend on the fixed 
value of Vx at which one measures the response. This is 
in contrast to the measurement of transverse response by 
trapping the polymer end at a constant Vy , while leaving 
it free to move in x-direction. Thus we reemphasize that 
the force-extension behavior depends on the kind of trap- 
ping potential used in an experiment. Apart from this, 
as we have shown, the orientational boundary conditions 
at the ends of a polymer and the ensemble of experiment 
will affect the force-extension behavior non-trivially. 
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Recently, using a Greens function calculation of a WLC 
polymer with one of its ends grafted, the presence of 
multiple maxima in p{vx,Vy) has been observedjs^. In 
this subsection, we have used MC simulations to study 
p{vx,Vy) for all the possible boundary conditions. We 
have shown that multiple maxima in p{vx,Vy) persists 
near t — A for all the three different boundary orienta- 
tions. We now utilize our theoretical methods as devel- 
oped in Sec|TT]to obtain the density plot for p{vx, Vy) at 
t = 4 for all three different end orientations (see FigfT4|). 
The results plotted in Fig[T4] using Nd — 5 already shows 
good agreement with Fig[T21 This clearly brings forth 
the presence of multimodality in the Helmholtz ensem- 
ble. Increase in the number of basis states Nd (infinite 
in principle) will lead to better agreement. 

Is it possible to test the results presented here experi- 
mentally? Fluoroscence microscopy of cortically confined 
actin filaments has been performed to extract their per- 
sistence length [ij. While most of the force extension mea- 
surements are usually done in 3D, however, experiments 
on cortical actin filaments are, we believe, possible and 
one may indeed test the predictions described above in 
such systems. In this context, it is instructive to note 
that, to achieve the parameter regime i ~ 4, for instance, 
in actin polymers that have persistence length 16.7/xm 
one requires contour length L ^ Qliim which is easily 
achievable experimentally (filaments as long as 100/im 
have been reported [l[). L can be changed chemically by 
the addition of enzymes. To measure the multimodality 
predicted in this paper, one can perform direct video mi- 
croscopy of the conformations of actin polymers confined 
in a cell of depth ~ l/im, practically restricting all fluc- 
tuations in third direction making the embedding space 
2D as in Ref . [l| . In such a setup one can also attach one 
of the ends of the actin molecules to one of the confining 
glass walls of the cell that contains them. Thus a setup 
as in Ref. [H may be used to obtain the probability distri- 
bution of end to end vector for a free polymer as well as 
a polymer with one end grafted. With actins of contour 
length ~ 67jum the probability distribution of end to end 
vector should show multimodality implying bistability in 
the force extension measurement in the Helmholtz en- 
semble. In this stiffness regime even the projected prob- 
ability distribution is expected to show multimodality in 
2D in contrast to 3D. In typical force-extension measure- 
ments one or both the ends of a polymer are attached to 
dielectric or magnetic beads to hold the ends optically or 
magnetically. In a recent study[l3| it was shown that to 
extract physically meaningful results from such experi- 
ments on dsDNA one has to incorporate bead geometry 
explicitly in the theoretical modelling, since the typical 
bead radius R is in between 0.05 - O.S/xmT^ which is 
about one to ten times the persistence length (A — 50nrn) 
of dsDNA. However, since the persistence length of actin 
filament is much larger (A = 16.7/xm), for actin R/\ is in 
between 0.003 - 0.03, one does not need to worry about 
bead geometry in analyzing the force-extension results for 
actin filaments. Our theoretical predictions can thus be 
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FIG. 15: (Color online) Representative conformations at 
t = A, for polymers with one end grafted in a;-direction, are 
shown. At this stiffness value the conformations are either 
localised corresponding to Gaussian chain like behaviour (a), 
or extended corresponding to rigid rod like behaviour (b). 



straightaway tested in experiments on actin. Note that, 
semifiexible polymers in 3D are also expected to show 
multimodality in the probability distribution of end to 
end vector and therefore non-monotonic force-extension 
in Helmholtz ensemble at least with end orientations free 
to rotateflBl- In Fig [15] we show some typical confor- 
mations of a semifiexible polymer at i = 4 lying in 2D 
embedding space. At this stiffness, there are two max- 
ima in the end to end separation, one is near the zero 
separation and the other is near the full extension (see 
FigsfT^ and [H)) . FiglT^Ka) shows some representative 
conformations with nearly zero extension and Fig|15l(b) 
shows the same at near full extension. 



V. CONCLUSION 

In this paper, we have shown that the results of force- 
extension experiments on semifiexible polymers would 
depend on ensemble, constraints on end orientations, di- 
mensionality of embedding space and the kind of trapping 
potential used. In an earlier work we have shown the 
presence of multiple maxima in the probability distribu- 
tion of end to end distance of a free polymer at interme- 
diate stiffnesses (near t ^ 4) that lead to non-monotonic 
force-extension in the Helmholtz ensemble [l6|. In this 
paper, we have demonstrated that though the details of 
the end to end distribution depends crucially on the con- 
straints imposed on the end orientations, the multimodal- 
ity in the distribution always persists. In this paper, we 
have used a mapping of the WLC model to a quantum 
particle on a sphere to obtain probability distribution 
of end to end separation and force-extension in various 
ensembles taking care of the particular types of the con- 
straints on the orientations of polymer ends. We have 
made a number of predictions about the end to end statis- 
tics and the force-extension behaviors in the Helmholtz 
and the Gibbs ensemble. We have used MC simulations 
against which we have tested the theoretical predictions 
and always obtained very good agreement. Experiments 
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using the laser trap to hold the ends of a polymer al- 
lows all possible end orientations, whereas the magnetic 
tweezers can be used to fix the end orientations and see 
the impact. On the other hand it is possible to obtain 
video microscopy of semiflexible polymers like actin to 
obtain the probability distribution of end to end separa- 
tion. Thus it is possible to test our theoretical predictions 
in experiments. In this work we have restricted ourselves 
to 2D. At the onset we have shown that an important 
feature of polymer statistics, multimodality in projected 
distributions, is dependent on the dimensionality of em- 
bedding space. In three dimensional free polymers, mul- 
timodality in projected probability distribution is impos- 
sible, however presence of this is a reality in 2D. We have 
shown that depending on whether the dimensionality d^, 
in which the trapping potential traps the polymer ends, is 
same or less than the dimensionality of embedding space 
d, the physically relevant Helmholtz free energy would be 
obtained from the probability distribution of end to end 
vector or a dr-dimensional projection ( [d — dr] dimen- 
sional integration ) of it. After projection, multimodality 
in the distribution function of end to end vector may or 
may not survive, thereby affecting the qualitative fea- 
tures of the force-extensions. Fixing the orientation of 
a WLC polymer at one end we have studied the pro- 
jected probability distributions in the longitudinal and 
transverse directions. The transverse fluctuations and 
the force-extensions found from our theory show excel- 
lent agreement with MC simulations in Ref.^3l- If ori- 
entations at both the ends are kept fixed, the polymer 
properties vary depending on the relative angle between 
the two grafted ends. For example, multimodality in pro- 
jected distribution depends on the relative angle. The 



full statistics of the WLC polymers are encoded in prob- 
ability distribution of end to end vector. Our simulations 
and theory have clearly shown that, the multiple maxima 
feature in this probability distribution in the intermedi- 
ate stiffness regime (near t — A) survives the fixing of 
end orientations. Similar studies in 3D remains to be an 
interesting direction forward. Multimodality in probabil- 
ity distribution may show multistability in the time-scale 
the end to end separation of a WLC polymer spends in 
each of the free energy minima. In polymer looping, the 
closing time and the opening time of the two ends of a 
free polymer depends on the polymer stiffness. The im- 
pact of the triple minima in the Helmholtz free energy on 
these time- scales of a free polymer remains to be stud- 
ied. This might be of importance in understanding the 
very fast time-scale of transcription, with respect to the 
diffusion time, in the process of gene expression [4|. We 
intend to report on some of these problems in future. 
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